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Synopsis 

This manuscript is intended to be a companion volume to the "METCAN User’s 
Manual" and the "METCAN Demonstration Manual". The primary purpose of the manual 
is to give details pertaining to micromechanics and macromechanics equations of high 
temperature metal matrix composites that are programmed in the METCAN computer code. 
The subroutines which contain the programmed equations are also mentioned in order to 
facilitate any future changes or modifications that the user may intend to incorporate in the 
code. Assumptions and derivations leading to the micromechanics equations are briefly 
mentioned. 


Introduction 

High temperature metal matrix composites offer great potential for use in advanced 
aerospace structural applications. The realization of this goal however, requires concurrent 
developments in (1) a technology base for fabricating high temperature metal matrix 
composite structural components, (2) experimental techniques for measuring thermal and 
mechanical characteristics, and (3) computational methods to predict their behavior. In fact, 
it might be argued that the development of computational methodologies should precede 
the others because the structural integrity and durability of high temperature metal matrix 
composites (HT-MMC) can be computationally simulated, and the potential payoff for a 
specific application can be assessed, at least qualitatively. In this way, it is possible to 
minimize the costly and time consuming experimental effort that would otherwise be 
required in the absence of a predictive capability. 

Recent research into computational methods for simulating the nonlinear behavior 
of high temperature metal matrix composites at NASA Lewis Research Center has led to 
the development of the METCAN (Metal Matrix Composite ANalyzer ) computer code 
(References 1-3). It is a computer code to simulate the high temperature nonlinear behavior 



of continuous fiber reinforced metal matrix composites. METCAN treats material 
nonlinearity at the constituent (fiber, matrix, and interphase) level, where the behavior of 
each constituent is modeled using a time-temperature-stress dependence. The composite 
properties are synthesized from the constituent instantaneous properties by making use of 
composite micromechanics and composite macromechanics models. Factors which affect the 
behavior of the composite properties include the fabrication process variables, the in-situ 
fiber and matrix properties, the bonding between the fiber and matrix, and/or the properties 
of the interphase between the fiber and matrix. The METCAN simulation is performed as 
point-wise analysis and produces composite properties which can be incorporated into a 
finite element code to perform a global structural analysis. After the global structural 
analysis is performed, METCAN decomposes the composite properties back into the 
localized response at the various levels of the simulation. At this point the constituent 
properties are updated and the next iteration in the analysis is initiated. This cyclic 
procedure is referred to as the integrated approach to metal matrix composite analysis and 
is depicted in Figure 1. 

METCAN can be used as a stand alone code to simulate the nonlinear behavior of 
MMC or can be used as a module in any structural analyses package to predict the metal 
matrix composite structural response. It can incorporate the fabrication induced residual 
stresses in its simulation of composite response. A variety of physical effects such as 
temperature dependence, stress and stress rate dependence, fatigue due to mechanical and 
thermal cyclic loads can be accounted for in METCAN analysis. Furthermore, the code can 
be updated with relative ease to include any other physical influences that one might want 
to consider in a specific analysis. The documentation for the code consists of three volumes: 
1) METCAN User’s Manual, 2) METCAN Demonstration Manual and 3) METCAN 
Theoretical manual. The details of how to use the code are discussed in the first volume. 
The various features of METCAN and the types of analysis that METCAN can peroform 
are illustrated through several demonstration problems in the second volume. The 
theoretical manual lists all the equations pertaining to the high temperature metal matrix 
composites with the idea of enabling the user to update or modify the source code if 
necessary. 


Scope 

The scope of the present manual is to summarize the micromechanics equations that 
are programmed in the METCAN computer code. The micromechanics equations provide 
a basis for expressing the equivalent ply level properties of the metal matrix composite in 
terms of its constituent in-situ properties. Three constituents are recognized in the present 
code. They are the fiber, matrix and the interphase. It should be noted that the interphase 
can also be replaced by a compliant layer. The only difference is by virtue of the definition 
the interphase is a region formed between the fiber and matrix as result of a chemical 
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reaction reducing the diameter of the fiber. A compliant layer, however, is a coating on the 
fiber provided intentionally to prevent such chemical reaction. Consequently, the effective 
fiber diameter is increased in the case of compliant layer. The mechanics involved in 
integrating the ply level properties to the laminate level properties (macromechanics) is well 
known and any text book on composites such as Reference [4] can be consulted. These 
details are not repeated here. 


Coordinate Axes 

METCAN adopts two sets of cartesian coordinate axes for presenting the results. They are 
the material axes and the structural axes as shown in Figure 2. The material axes are defined 
with respect to fiber orientation. The direction along the fiber is denoted as 1-1 or 11 and 
the directions transverse to the fiber are denoted 2-2 (22) and 3-3 (33) respectively. All the 
ply level properties and responses are expressed with respect to the material axes system. 
The global/laminate level properties and responses are given in the structural axes system. 
They are denoted by x, y, and z. Any quantity with a normal orientation is represented by 
the subscripts xx, yy, and zz. The shear orientations are indicated by subscripts xy, yz, and 
zx.. Also, the coordinate axes definitions are printed as a part of the output for easier 
interpretation. 


Units 

The primary units chosen in METCAN are the British system of units. All of the constituent 
properties in the resident data bank are expressed accordingly: inches (for length), pounds 
(for forces), degrees Fahrenheit (for temperatures), BTU s for heat, seconds/hours (for time) 
and psi/ksi/Mpsi (for stress and stiffness related quantities). A table of units (Table 1) is 
provided (by default as one of the output options) at the beginning of the METCAN output 
for easy reference. Any user wishing to employ another system of units may do so by 
updating the data bank properties into the preferred system. The output will then reflect 
the appropriate unit system chosen by the user. Nevertheless, the user should exercise 
some caution as there may have been some hard wired conversions of units within the 
source code (e. g. conversions of psi to ksi, or mpsi are done by dividing with 1000 or 
1000000 with in the code in order to print some of the variables in the most appropriate and 
concise format) 
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Glossary of symbols 


11, 22, 33 

subscripts for longitudinal and transverse 
directions in the material coordinate 
system 

f, m, d, and 1 

subscripts for fiber, matrix, interphase 
and ply 

A, B, C 

unit cell subregions 

C 

heat capacity 

E 

Young’s modulus 

D 

Current diameter of the fiber 

D c 

Outer diameter of the fiber 

G 

Shear modulus 

K 

thermal conductivity 

T 

temperature 

0 

suffix for the reference (room temp. ) 
quantities 

k 

fiber volume ratio 

a 

thermal expansion coefficient 

V 

Poisson’s ratio 

a 

stress 

P 

density 

AT 

temperature differential 


Nonlinear Material Behavior Models 

METCAN treats material nonlinearity at the constituent level through the incorporation of 
a multifactor interaction relationship (MFIR). The MFIR models the material behavior in 
the form of a time-temperature-stress dependence of the constituent properties. The MFIR 
and the rationale for its utilization are expressed in Figure 3. The form of the MFIR is a 
product function of terms raised to specified powers that define the values of the material 
property (P) of the fiber, matrix, and interphase. The constituent properties (P) which are 
determined through the MFIR are modulus, Poisson’s ratio, strength, coefficient of thermal 
expansion, thermal conductivities. The density and heat capacity are considered as constants. 
For each fiber, matrix and interphase property, there is generally a distinct set of values of 
the exponents necessary to express the functional definition of the property. The various 
exponents are determined from the appropriate experimental data. The variation of a 
typical property as a function of one such factor is shown in Figure 4 for various exponent 
values. As seen from the figure, one can nullify any of the physical effects (e. g. temperature 
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dependence) by simply choosing a "zero" exponent. A Unear variation can be achieved by 
taking a value of "one" for the exponent. The property value must be matched carefully at 
the reference point and at the final point from the appropriate experimental data. The 
actual nonlinear behavior between these two points is dependent on the value of the 
exponent which is usually of a second order effect. 

The coding for the material behavior modeling is done through a set of six subroutines in 
the METCAN code. The thermal properties are updated in the subroutines THERF, 
THERM, and THERD respectively for fiber, matrix and interphase properties. The 
mechanical properties are updated in the subroutines MECF, MECM, and MECD 
respectively for fiber, matrix and interphase properties. 


METCAN unit cell 

In METCAN the global behavior is obtained by successive integration of the 
response/properties from the constituent level to laminate level using micromechanics and 
macromechanics. The micromechanics embedded in METCAN is based upon the usual 
assumptions pertaining to the mechanics of materials approach [Ref. 4-5]. The derivation 
of equivalent composite properties using micromechanics theory starts with the identification 
of a representative volume element or unit cell. It is the smallest region or piece of material 
over which the stresses and strains are macroscopically uniform. However, within the unit 
cell the stresses and strains are nonuniform due to the heterogeneity of the material. It is 
assumed that these cells are arranged in a regular square array pattern. The principal 
assumptions involved in the mechanics of materials approach are (1) fiber and matrix are 
subjected to the same strain in the fiber direction of a unidirectional fibrous composite and 
(2) the same transverse stress is applied to both the fiber and matrix in the direction 
transverse to the fiber. The details of the representative volume element chosen in the 
development METCAN is shown in Figure 5. The unit cell for METCAN consists of a 
maximum of three regions, A, B, and C. The region A consists of matrix only, the region 
B consists of a combination of matrix and interphase, and the region C consists of a 
combination of all three constituents: fiber, matrix and interphase. Note, that in the 
following discussion the interphase is treated as a separate constituent with distinct 
properties. Thus, it can either represent a zone formed due to a chemical reaction between 
the fiber and matrix or a separate layer provided intentionally to prevent such a reaction. 
In the figure 5 D represents the current diameter of the fiber and D 0 represents the outer 
diameter. The outer diameter is either the undamaged fiber diameter if there is an 
interphase or the diameter of the fiber with the compliant layer coating if there is a 
compliant layer. The different regions A B, and C of unit cell facilitate the representation 
of nonuniformity in the local stress distribution. 
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Microstresses in Regions A, B, and C of METCAN unit cell 
The properties of the constituents forming the unit cell in the regions A, B, and C are 
updated incrementally at the end of each iteration depending upon the local response 
(stress, stress rate, temperature etc..). The local response is obtained by progressively 
decomposing the global response to ply response (ply stresses) and applying the 
micromechanics at the unit cell. The stresses in the various regions of unit cell are termed 
as "microstresses" in METCAN. These calculations are performed in the METCAN 
microstress calculation module which consists of a set of four subroutines STRESF, 
STRESM, STRESD and STRESI. The first three pertain to the calculation of microstresses 
in the various regions for fiber, matrix and interphase respectively. Operations common to 
these three are grouped together into a common subroutine STRESI which is called before 
calling the rest of the three subroutines. The following few paragraphs lists the various 
equations that are currently programmed in METCAN to perform the microstress 
calculations. 


Stresses due to temperature differences 

The following equations show how microstresses due to temperature differences in the fiber, 
microstresses in matrix in regions A, B and C, and microstresses in interphase in regions B 
and C are calculated in METCAN. It should be noted that temperature differences cause 
microstressess due to two reasons. The first one is due to differences in constituents thermal 
expansion coefficients. The second one is due to the ply stresses that are developed as a 
result of different ply orientations. The latter is described in the next section "stresses due 
to mechanical loads". 

Microstresses along the longitudinal (11) direction: 

Fiber Microstresses: The fiber longitudinal microstress due to a uniform temperature 
difference AT can be given by 

°m = {«iu _a fii} E m AT (1) 

In the above equation the longitudinal stress in fiber due to a temperature difference is 
calculated using the differential expansion between the fiber and ply. However, in 
METCAN it is calculated from the stress equilibrium along the fiber direction using the 
following equation. 



The above equation forces fiber stress to be in equilibrium with those of matrix and 
interphase. 

Matrix Microstresses: The matrix microstress has two components a mllT and a mllNU where 
the part of the subscript T and NU denote the normal and Poisson’s contributions 


6 


respectively. They are given by 

°mllT = { a iU -a inll^mn^T 

°mllNU = [ ^ { V m21A( a t22 _a m22) +V m21B( a i22 _<X m22) + 
V m 2 ic( a l 22 _ a n 22 ) + V m3XA( a l22 ” a m33) + 
V m31B( a i22 _a xii33) +V iii31c( a l22 -a in33) } 


( 3 ) 

( 4 ) 


°mll = °mllT + °mllNU 


( 5 ) 


Interphase Microstresses: The interphase microstress, similarly, has two components a dUX 
3. rid ^diiNu are given by 

°dllT = { a lU -a dll}®dll^^ ^ ^ 


J dllNU 


= f l/4{v d21B (a 122 _ ad22) +V d21c( a i22 _a d22) + 

V d31B( a l22 -a d33) +V d31A( a l22 _a <>33) } pdll A '*' 


'dll 


dUT + °dl!NU 


Microstresses along the transverse (22 and 33) directionsi 
Transverse fiber microstresses: 


°f22C a m22C 
°f33C = °m33C 


( 7 ) 


( 8 ) 


( 9 ) 


Note that in the above it is assumed that both the fiber and matrix develop the same 
mi crostress in 22 and 33 directions. 

Transverse (22 &33) thermal microstresses in matrix Region A: The tranverse thermal 
microstresses due to a temperature differential AT is computed by considering the 
differential expansion of the region A and the ply. The end result is 


°m22A = { a l22 _a m22A +V ml2A( a Lll _a mll) +V m32A( a B3 a in33A)} ATE in22A 


( 10 ) 


°m33A = { a l33 _a m33A +V iiil3A( a lll _a mll) +V in23A( a l22 a m 22A)) ATE m33A 


(ID 


Transverse (22 & 33) thermal microstresses in matrix Region B: These microstresses are 
derived along the same lines as the above. However, since the region B consists of both 
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matrix and interphase, there are contributions from both constituents. Furthermore, both 
interphase and matrix will have the same stress given by the following equations: 


J tn22R 


i/5 

a -ft,) 

d-^) 

C C2 


V <U2B^ a Ul -a dll^ “ V ml2B^ 


V d32B^133~^d33B^ _ V m32B^ ®133 ®m33B^ 


m33B 


(1-^) (1-^) 




d22B m22B 


ATfl-^E^ 




-<t22B 


"m22B I 


a 


(l-V^) 
133 


V dl3Bt a Ul“ a <»nJ “ V ml3Al a Ul -a mlll 


V d23B^ g m~ g Hy»p} V m23B^ ^122 _ ^m22B^ 


v*r 


(1-^) (1-^) 


a d33B a m33B 


AT(1 y/^f)E (1 33 B 




^<I33B 


"m33B 


( 12 ) 


(13) 


Transverse (22 & 33) thermal microstress in matrix region C: These are calculated by 
forcing the equilibrium in the respective directions. The following are the appropriate 
equations: 


°m22C = 
°m33C 



(14) 

(15) 


Transverse microstressess in the interphase regions B and C: These are assumed to be same 
as those for the matrix regions B and C. 


°d22B = °m22B 
°d33B = a m33B 
°<J22C = °m22C 
°d33C = °m33C 


(16) 


This completes the summary of equations for microstresses due to temperature differential. 
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Stresses due to mechanical loads 

The mechanical loads applied globally are progressively decomposed using 
macromechanics/laminate theory to ply stresses and these stresses are passed on to the 
microstresses calculation module. It should be noted that, eventhough they are treated as 
mechanical loads, contributions arising due to thermal loads are included in the ply stresses. 
The ply stresses are assumed to have been applied to the unit cell. The constituent 
microstress contribution corresponding to each of the these ply stresses in each of the 
regions A, B and C are calculated using the mechanics of materials approach. 


Fiber Microstresses due to applied plv norm al stresses 

Stresses along the longitudinal (11) direction: The fiber longitudinal (11) stresses due to the 
three normal applied ply stresses are given by 

o, 


(1) _ n 

°m " fc ni 


■*111 


■"in 


( 2 ) 

°ni 


„0) 


■mi 


/ \ ^fii 

" ( v l21 _V Cl)'^ °122 

*"122 


(17) 


~ ( v m v Gi) 


"m 


01 G1 J 

off, * ogl 


■’133 


"03 

+ °fll 


111 

Note that in the above the superscripts and ^ represents the contributions to a 

specific microstress due to the applied ply stress in 11, 22 and 33 directions respectively. 
Stresses along the transverse (22) direction: The fiber transverse (22) microstress due to the 
three normal applied ply stresses are given by 


°f22 - E22C p 

^122 


122 



V 112 V ml2C - ^f 



dl2C 





(18) 

(19) 

( 20 ) 


In the above the E^c is given by 
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E22C - 


0-/5 t 

f L°1 

+ &R 

Em22C 

Ed22c[ D 0 J 

Ef22 D 0 


( 21 ) 


J C2 


- 4 . . ^(3) 

“ Of2 2 + ^£22 + °f22 


( 22 ) 


Stresses along the transverse (33) direction: These are analogous to the 2-2 direction stresses 
and are given by 


_(3) _ TJ IM 
°f33 - =330 VT~ 
M33 


133 


(23) 


„<!) - 
°£33 " 


'U3 




'ml3C { 1_ D 


■ V «U3C + 7T V fl3 

U 0 / j 


u lll 

Bui 


E33C 


(24) 


( 2 ) _ _ 
°Q3 - 


D f 


? 123 (l - ^) V m23C~^f I j * 


' V <C3C + _ pr V «22 
U 0 


P 

P “330 
*122 


(25) 


In the above the is given by 

E33C = 


(1 -fe f) ^ yjkf 


. A D 

^m33C ^d33C 

l D «l 

8 W 

PI 


(26) 


°f33 ~ °f33 + °f33 + °I33 


(27) 


Matrix microstresses due to applied plv normal stresses 

Stresses along the longitudinal (11) direction in region A: The matrix longitudinal stresses 
due to the three normal applied ply stresses are given by 


m _ t, U ui 
°mllA ” “mllA'i 


L Em 


«&. ■ (v E1 -v G ,)^i<. *• 


ml 1A V *121 p 122 ^ 


(28) 


E m kf 

°milA = ( V 131 _V f31^B °B3: 


. 0 ) 


£133 “(1-ty 
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Stresses along the longitudinal (11) direction in regions B and C: These are given by similar 
expressions as above. However, the Poisson’s contributions (the second and the third of the 
above equations 28) are taken as same as above. 

a, 


_(i) _ p 

°mllB “ ^mllBrj 

“U1 


111 


( 29 ) 


( 1 ) _ £ 

°ml 1C ml 1C tj 

nil 


111 


The total matrix microstress is given by adding each contribution. 


°mllA 0 


°mlic 0 


(1) 


(2) 


n O) 

mllA 

+ 

°mllA 

+ 

°mllA 

(1) 


(2) 


_(3) 

f mllB 

+ 

OmllB 

+ 

a mllB 

(1) 


Q) 


_(3) 

r mllC 

+ 

°mllC 

+ 

°mllC 


( 30 ) 


Stresses along the transverse (22) direction in regions A, B, and C: The transverse 
microstresses in the matrix region B and C are calculated first and the region A microstress 
is obtained by forcing equilibrium in 22 direction. The microstresses in region B due to the 
three normal applied ply stresses are given by 

(2) _ ^22B 

°m22B = a i22“ 


AD 

°m22B 


^122 
^11 

3^28 = " ‘^ E ^ {v i32 -V m32B (l -^-^Vds^} 

^12 

^22B = 


( 31 ) 


D d22B 


(fit + (1-v^) 


~'d22B 


~'m22B 


) 


By adding the above three contributions the total microstress in region B is obtained as 

_ J» . „» * ( 32 ) 


°m22B “ °m22B + ° m22B + °m22B 

The region C stresses in matrix, similarly, can be expressed by the following set of equations 
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where is given by equation (21). 

By adding the above three contributions the total microstress in region B is obtained as 


_ _ m . a) . o) 

°m22C " °m22C + °m22C + °m22C 



Stresses along the 33 direction in Regions A, B, and C: The microstresses in 3-3 direction 
are given by analogous expressions as those for the 2-2 directions and may be obtained by 
appropriately changing the subscripts 22 to 33. For the sake of completeness they are given 
below: 
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_(2) _ _ ^33B 

°m33B °133“ 


^133 


a (1) 

a m33B 


= p E33B { V tl3- V ml3B(l ~\/^f)~V^f V dl3B} 
^lll 

®m33B = ^^ ~ ~ { V 123 ~ V m23b( 1 ~^ V <C3b} 


( 36 ) 


^33B 


”'d33B 


(fet + (1 \^f) 


-d33B 


-'m33B 


By adding the above three contributions the total microstress in region B is obtained as 


.(i) 


ri) . „(3) 


°m33B " °m33B + °m33B + °m33B 


(37) 


The region C stresses in matrix, similarly, can be expressed by the following set of equations 


(2) _ E 33C 

°m33C " °133 t? 

n33 


d33C * ( V «* 

l _V m23B^ ( V 423C | 1 "'5^| + '5^' 


(2) _ 0 122^33C 

°m33C » 


fl3 


£23 


(38) 


where is given by equation (26). 

By adding the above three contributions the total microstress in region B is obtained as 

(1) _(2) (3) (39) 

°m33C = °m33C + °m33C + °m33C 


The Region A stresses are given by 
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(3) _ J 

a m33A ” °133" 


I (1) 

J m33A 


- / m33A 

Ej33 


J m33A 


_ _ _G) , „(?) 

°m33A “ °m33A T °m33A 


(1) . 

D (1) 

°m33B + 

— a m33C 

u o 1 

(2) 

D (2) 

°m33B + 

■ZT" G m33C 

D 0 

4. „0) 


+ °m33A 



( 40 ) 


Interphase microstresses due to applied plv normal stresses 

Stresses along the longitudinal (11) direction in Region B: The interphase longitudinal 
stresses due to the three normal applied ply stresses in the Region B are given by 

o, 


_0) _ p 

°dllB ~ “dllB^ 


'111 


"Ul 


°!uiB = " ( V 121 “ V «121 b) 


"dllB , 


"122 


J 122 


(41) 


°dllB ( v 131 V d 31 B)“S °133 

“133 

_ - -CO 4. „P) * „0) 

°dUB “ °dllB + a dllB + a dUB 

Stresses along the longitudinal (11) direction in Region C: These are obtained by changing 
the subscript B to C in the above equations. Accordingly, these are given by 

o, 


(1) _ tj 

°dllC “due" 


'U1 


-'ll! 


„(2) _ \ d!lC_ 

°dlic ( v m V d21C' - p °I22 

“l22 

°dllC = ~ ( V 131 ~ V d31C^ - p °133 

“l33 


(42) 


= -CO 


.(2) 




°dllC “ a dllC + °dllC + °dlIC 


Stresses along the transverse (22 and 33) direction in Regions B and C: These are same as 
the matrix microstresses in the regions B and C. 
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( 43 ) 


®d22B ” ®m22B 
°d22C = °m22C 
°d33B = °m33B 
°d33C = °m33C 

The above completes the discussion on constituent microstresses in various Regions A, B, 
and C It should be noted that the mircrostresses due to temperature differential must be 
added to those due to mechanical loads to obtain the combined microstresses due to 
thermo-mechanical loading. 

Apart from these microstresses which are basically the stresses due to normal loads, there 
exist the shear stresses due to the applied shear force resultants. These are also obtained 
in a s imil ar fashion and are given in the following paragraphs. 


Fiber microstresses due to applied ply shear stresses 

The fiber microstresses due to the three applied ply shear stresses are given by 


'112 j 


J fl2 


'112 


'fl2 


'ml2 


'fl3 


[ 

D 0 


„ 

®113 

G U3 

b 

_ D G ml3 
D 0 G dl3 


G f23 

°123 r . 

G 123 


'fl2 






J f23 




(44) 


Matrix microstresses due to applied plv shear stresses 
Region A: The matrix microstresses due to the three 
by 


applied ply shear stresses are given 


°ml2A 


°112 


G ml2A 

°112 


0 ml3A 


'113 


'ml3A 


'113 


» °m23A 


°123 


G m23A 

G 123 


(45) 
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Region B: The microstresses in region B due to the three applied ply shear stresses are 
calculated by 


°ml2B " °112' 


*ml2B 


■*112 




°ml3B °113' 


J ml3B 


*113 


1 




(46) 


°m23B ~ °123" 


*m23B 


*123 


Region C: The region C stresses are calculated by 



Interphase microstresses due to applied ply shear stresses 
Region B: The region B interphase stresses are calculated from 


(47) 
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°dl2B 


'112 


^d!2B 

(■*112 



°dl3B 



°d23B 



( 48 ) 


Region C: The region C stresses are given by 


°dl2C °U2 r 


'dl2C 


'112 


i K~ ( t rt ® \ ^b 12 D ^ml2^ 
V ^0 U dl2 ^0 U fl2 ) 


a dl3C °U3 


'dl3C 


■*113 


\ 1 


i FT (i /i ^ml3 D ^ml3 


°d23C “ °123 r 


J d23C 


J t23 




( 49 ) 


This completes the equations for microstresses due to applied ply shear stresses. Overall, 
there are six fiber microstresses (three normal and three shear), 18 matrix microstresses 
(nine normal and nine shear), and 12 interphase microstresses (six normal and six shear) 
that result due to the six ply stresses (three normal and three shear). The ply stresses are 
calculated in the subroutine COMSA using classical lamination theory for each incremental 
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load step. These are then passed to the microstress module comprised of the subroutines 
STRESI, STRESF, STRESM and STRESD for the computation of the incremental 
microstresses which are cumulatively added to obtain the total microstress in each 
constituent and in each region. These microstresses are input to the MFIR module for 
calculating the stress dependence of the constituent material behavior. Also, at the end of 
each computation, these microstresses are compared with the respective allowable strengths 
(e. g- CT ni k compared with S fuxc ) to check for any failures in the constituents. If a 
particular constituent shows stress failure, its respective stiffness is equated to a negligible 
value (e.g. E ni = 0.01 psi) in order redistribute the load and change its path. 

Micromechanics 

Composite micromechanics theory refers to the collection of physical principles, 
mathematical models, assumptions and approximations employed to relate the behavior of 
a simple composite unit (e.g., lamina or ply) to the behavior of its individual constituents. 
The primary objective of composite micromechanics is to determine the equivalent 
mechanical, and thermal properties of a composite ply in terms of the properties of the 
constituent materials. The mechanical properties refer to normal moduli, shear moduli, and 
Poisson’s ratios. The thermal properties include the thermal conductivities, heat capacity and 
the coefficients of thermal expansion. The derivation of equivalent composite ply properties 
using micromechanics theory starts with the identification of a representative volume 
element or unit cell as shown in Figure 5. It is the smallest region or piece of material over 
which the stresses and strains are macroscopically uniform. However, within the unit cell the 
stresses and strains are nonuniform due to the heterogeneity of the material. The unit cell 
can consist of a fiber, matrix and/or an interphase. It is assumed that these cells are arranged 
in a regular square array pattern. Equivalent properties for the ply are then derived in terms 
of the constituent material properties based on the mechanics of materials approach. Other 
assumptions involved in this approach are: 1) fiber and matrix are subjected to the same 
strain in the fiber direction of a unidirectional fibrous composite, and 2) the same transverse 
stress is applied to both the fiber and matrix in the direction transverse to the fiber. 

Physical Properties 

The mass density of the ply is assumed to be given by the rule of mixtures type equation as 
shown below: 

Pi = fcfPf + ^nPm + (i-VkJPd < 50 ) 

The above is calculated in the subroutine PLYMAT. 

Mechanical Properties 

The micromechanics equations which facilitate the computation of ply equivalent properties 
are programmed in the subroutine MECHPL in METCAN. The actual details of derivations 
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for the ply longitudinal and transverse moduli are given in the Appendix. The derivation of 
the remaining properties follow along the same lines as described in the Appendix. The 
following is a brief list of the various equations: 


Ply Normal Moduli: The longitudinal modulus E, u of the ply is calculated based on the 
assumption that the compatibility of longitudinal displacement requires equal strains for the 
composite and constituents [1]: 


Bui ~ + kf 



(50) 


The transverse normal moduli and are derived by first computing the equivalent 
moduli in each of the regions A, B and C and then combining them by assuming that the 
subregions act as parallel elements when subjected to a transverse load. To compute the 
equivalent moduli of each subregion it is assumed that the stress due to transverse load is 
same in each of the constituents. This leads to the following expressions for the transverse 
moduli: 



K' 




(51) 


(52) 


Ply Shear Moduli: 

The derivation for the shear moduli in 1-2 and 1-3 directions follow along the same lines as 
that of the transverse modulus. The shear modulus in 1-2 direction is given by 



The shear modulus in 1-3 direction is given by 
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(54) 
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The shear modulus in 2-3 direction is calculated by utilizing the transverse isotropy in the 
2-3 direction: 


G 123 
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(55) 


Alternatively, shear modulus in the 2-3 direction can also be derived independently using the 
same criteria as that of transverse modulus derivation. It is given by 
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Note that currently in METCAN this equation is not used. 


Poisson’s Ratios: The derivation for the major Poisson’s ratio in 1-2 and 1-3 direction follow 
the derivation for the longitudinal modulus where in the fiber and matrix elements are 
assumed to be acting in parallel. The respective equations are 
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(58) 


Poisson’s Ratio in 2-3 direction: The Poisson’s ratio in 2-3 direction is derived by using a 
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combination of series and parallel elements of the various subregions A, B and C. The 
equation is given by 
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Thermal Properties: 


The properties under this category are the three thermal expansion coefficients, the three 
thermal conductivities and the heat capacity. These are computed in the subroutine 
THERMPL in METCAN. 

Longitudinal thermal expansion coefficient: This is derived based on the observation that 
free expansion in 1-1 direction produces no resultant load in that direction. This implies that 
the integrated effect over the cross section of the sum of the individual stresses developed 
in each constituent due to the mismatch of the thermal expansion coefficients must add up 
to zero identically. This leads to the following expression for the longitudinal thermal 
expansion coefficient: 


a lll " a mll + M 
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Transverse thermal expansion coefficients: The derivation for the transverse thermal 
expansion coefficients is carried in two stages analogous to the derivation of the transverse 
moduli. First, equivalent thermal expansion coefficients are written for each subregion A, 
B, and C. During this stage the individual constituents are assumed to act in series. The 
overall equivalent thermal expansion for the unit cell is then written by assuming the 
individual regions to act in parallel. The following are the equations for the two transverse 
thermal expansion coefficients: 
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Longitudinal heat conductivity: The longitudinal heat conductivity of the ply is derived by 
assuming that the three constituents are in parallel which yields the following rule of 
mixtures type expression: 

K U1 - (<3) 

Transverse heat conductivities: To derive the transverse conductivities the approach adopted 
for the derivation of the transverse moduli is followed. The equivalent heat conductivity for 
each subregion A, B and C is first expressed assuming that the constituents are acting in 
series followed by the expression of the ply conductivity assuming that each subregion is 
acting in parallel. This yields the following expressions for the transverse conductivities: . 
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Heat capacity: The equivalent ply heat capacity is obtained by the simple rule of mixtures 
type equation. 
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where p, is the equivalent ply density. 


Macromechanics 

Macromechanics is the study of composite behavior wherein the material is presumed 
homogeneous and the effects of the constituent materials are detected only as averaged 
apparent properties of the composite. Thus, a collection of procedures which allows to 
proceed from the lamina or ply level properties to the global or laminate level properties 
where the heterogeneities associated with fiber, matrix and interphase are completely 
smeared out constitutes the theory of macromechanics. These procedures utilize the 
classical laminate theory to integrate the ply level properties and obtain the laminate level 
properties. These procedures are well known and are explained thoroughly in many text 
books. Therefore they are not repeated here. References [4-5] may be consulted for the 
theoretical details. Also, Reference [6] may be consulted for how these procedures are 
incorporated in METCAN, as many of these subroutines are taken from the earlier 
computer code ICAN which was primarily developed for polymer matrix composites. The 
principal subroutines that perform the macromechanical integration in METCAN are 
GACD3, GPCFD2 and COMSA. These are essentially the same subroutines that are used 

in ICAN [6,7]. 
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Incremental Iterative Solution Scheme 

The nonlinear thermo-mechanical response of high temperature metal matrix composites is 
obtained using an incremental iterative solution scheme as shown in the flow chart in Figure 
6. During each load increment a linear laminate analysis is performed. At the end of each 
such analysis the response is saved and compared with that of the previous step till 
convergence to a prescribed tolerance level is achieved. The properties of the constituents 
are varied in each iteration based upon the current response and temperature using the 
multi-factor-interaction-relationship. A different flow chart is shown in figure 7 which 
indicates the principal subroutines of the program that are active in each of the blocks 
shown in figure 6. Note that in this flow chart the modules performing output operations 
are intentionally removed to avoid clutter. Only the basic operations pertaining to nonlinear 
material model, micromechanics and macromechanics are indicated. This chart should be 
consulted along with figure 2 by those wishing to debug/update/modify the code. 

Convergence Criteria 

The convergence criteria used in METCAN is based on the global/laminate level strains. 
During each incremental load step analysis, the global incremental strains are preserved and 
are compared between successive iterations. If the differences are with in five percent 
convergence is assumed to have been achieved. The tolerance level is hard wired in the 
code. However, it can be easily altered. Also, alternative convergence criteria based on ply 
level strains/stresses/properties can be set with out much difficulty. 

Concluding Remarks 

The computer code METCAN is developed for simulating the nonlinear behavior of high 
temperature metal matrix composites. The code documentation consists of a set of three 
manuals. They deal with the code usage, demonstration of the code capabilities and the 
theoretical aspects pertaining to the high temperature metal matrix micromechanics. The 
present manual provides a summary of all the theoretical details embedded in the computer 
code METCAN. The emphasis, however, is placed specifically on the micromechanics, 
microstresses and the nonlinear material behavior models that are incorporated in 
METCAN. The primary intention is to provide information to the code user who might 
want to modify, update, and debug the code. Assumptions leading to the derivation of 
various equations are briefly mentioned with no details on the actual procedure. The details 
of the derivation is beyond the scope of the present manual. The actual usage of the 
computer code is discussed in the companion volume "METCAN User’s Manual. The 
various features of the code and how they can be utilized are discussed in the companion 
volume "METCAN Demonstration Manual". For a complete understanding of the code and 
its capabilities the user is advised to refer to all the three volumes. 
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Appendix 


In order to demonstrate the formal procedure involved in the application of composite 
micromechanics theory, derivations of the equations for ply normal moduli (E, n and E^) 
are explicitly developed below. The particular approach taken here relies on the principles 
of force equilibrium and displacement compatibility as defined from elementary mechanics 
of materials approach. 


longitudinal Modulus . 

Consider a square array unit cell model (see Fig. 5) subjected to a uniaxial load in the 
longitudinal direction. The equivalent composite (ply) load is defined from force equilibrium 
to be the sum of the constituent loads as follows: 

p t = p, + p„ ♦ p. (1) 


In the integrated average sense, Eq. (1) is rewritten as 

OjAj = a t Aj + o d A d + o m A m ( 2 ) 

where A represents cross sectional area. Division of the above equation by A, followed by 
the substitution of the respective ratios of the areas with the volume fractions of the 
constituents the following can be written: 

0! = Ofk/ + O d k d + ©An (3) 

where k/ is the actual fiber volume fraction based on the current fiber diameter D. 

Since the compatibility of longitudinal displacement requires equal strains for the composite 
and constituents («, = e f = e d = € m ), Eq. (3) can be differentiated with respect to strain 

to give 




do 


dej 


- f k 




(4) 


The quantities (d 0 /de) represent the slopes of the corresponding stress strain curves for the 
composite and constituents and in this context define instantaneous of tangent moduli. 
Hence, Eq. (4) becomes 

E U1 - E*' * E d k d ♦ E„k„ (5) 

The volume fractions can be expressed in terms of the inner and outer diameters of the fiber 
D and D c and the apparent fiber volume ratio k f .. Note that the apparent fiber volume 
ration is based on the original/outer diameter of the fiber D e . 
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( 6 ) 



With the aid of equations (5) and (6) the following expression can be developed for the ply 
longitudinal modulus: 








(7) 


Transverse Normal Modulus in 2-2 direction 


Consider the square array unit cell model again except that the fiber and interphase are of 
equivalent square cross section such that linear dimensions can be defined as follows: 


and 



( 8 ) 


s f = % ; s d = a^-aj ; ^ = a! (9) 

Assume a uniaxial load in the transverse direction and neglect Poisson"s effects. For the 
subregion C displacement compatibility yields 

s i € i = + S d € d + S m € m ( 10 ) 

and the force equilibrium results in equal stresses for the composite and constituents (a, = 
a f = a d = a m ). Hence, Eq. (10) can be differentiated with respect to stress to give 



(ID 


The quantities (de/d a) represent reciprocals of the slopes of the corresponding stress-strain 
curves for the composite and constituents and in the same context as before define 
reciprocals of instantaneous or "tangent" moduli. The above equation can be rearranged to 
yield the following: 
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( 12 ) 






+ 




With the aid of equations (7), (8) and (12) the equivalent modulus for the subregion C can 
be shown to be given by 



The equivalent modulus for subregion B can be deduced from Eq. (13) by letting D/D 0 
equal to unity. The result is 



(14) 


The equivalent modulus for subregion A is simply the matrix modulus or 



( 15 ) 


The ply transverse modulus E 122 is defined by assuming that subregions A, B, and C act as 
parallel elements when subjected to a transverse load. This is analogous to the case for E, n 
where the constituents are assumed to act in parallel. Hence, from Eq. 5 it can be deduced 
that 


Eft = E 1 C s f + E 1 \ + E^ (l 6 ) 

Division of the Eq. 16 by s, , followed by the substitution of the equations (7) and (8) and 
the results from equations (13 - 16) leads to the following final expression for the transverse 
ply modulus E^. 



Note that in the above the subscripts f, m and d are replaced by ”f22", "m22" and "d22" to 
indicate the direction of the loading. 
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Table 1. METCAN units for constituent, ply and laminate properties 


Property 

Symbol 

Unit 

ELASTIC MODULUS 

E 

psi 

SHEAR MODULUS 

G 

psi 

POISSOMS RATIO 

V 

non- dim 

THERM. EXP. COEPP. 

a 

in/in/P 

DENSITY 

p 

lb/in* *3 

FIBER DIAMETER 

D 

in 

HEAT CAPACITY 

C 

BTU/lb/P 

HEAT CONDUCTIVITY 

K 

BTU-in/HR/in**2/F 

STRENGTH 

S 

psi 

THICKNESS 

T 

in 

DISTANCE TO MIDPLANE 

Z 

in 

ANGLE TO AXES 

TH 

degrees 

TEMPERATURE 

T 

P 

STRAIN 

€ 

in/in 

STRESS 

a 

psi 

MEMBRANE LOADS 

N 

lb/in 

BENDING LOADS 

M 

lb-in/in 

FIBER VOLUME RATIO 

kf 

non- dim 
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Figure 1 .—Integrated multi-scale approach to metal matrix composite analysis. 
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(a) Single fiber. 



(b) Ply material coordinate system. 
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deg 
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(c) Laminate structural coordinate system. 

Figure 2.— Different coordinate systems used in 
METCAN. 
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Rationale: 
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• Gradual effects during most range, rapidly degrading near final stages 

• Representative of the in situ behavior for fiber, matrix, interphase, coating 

• introduction of primitive variables (PV) 

• Consistent in situ representation of alt constituent properties in terms of PV 

• Room-temperature values for reference properties 

• Continuous interphase growth 

• Simultaneous interaction of all primitive variables 

• Adaptability to new materials 

• Amenable to verification inclusive of all properties 

• Readily adaptable to incremental computational simulation 

Notations: 

P - property; T - temperature; S - strength; R - metallurgical reaction; N - number of cycles: 
t - time; over dot - rate; subscripts: O - reference; F - final; M - mechanical; T - thermal 

Figure 3. —Assumed multifactor interaction relationship to represent the various factors which influence in 
situ constituent materials behavior. 
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Figure 4. — Explanation of the multifactor interaction relationship. 


Figure 5. — M ETC AN unit cell showing the different regions A, B, 
and C of heterogenity and the idealization for the derivation of 
the micromechanical equations. 
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Figure 6. — M ETC AN flow chart 
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Figure 7. — METCAN flow chart showing the principal micromechanics and macromechanics routines. 
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